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We present a IIB supergravity solution dual to a spatially anisotropic finite-temperature N = 4 
super Yang-Mills plasma. The solution is static and completely regular. The full geometry can be 
viewed as a renormalization group flow from an ultraviolet AdS geometry to an infrared Lifshitz-like 
geometry. The anisotropy can be equivalently understood as resulting from a position-dependent 
#-term or from a non-zero number density of dissolved D7-branes. The holographic stress tensor 
is conserved and anisotropic. The presence of a conformal anomaly plays an important role in the 
thermodynamics. The phase diagram exhibits homogeneous and inhomogeneous (i.e. mixed) phases. 
In some regions the homogeneous phase displays instabilities reminiscent of those of weakly coupled 
plasmas. We comment on similarities with QCD at finite baryon density and with the phenomenon 
of cavitation. 



1. Introduction. The realization that the quark-gluon 
plasma (QGP) produced in heavy ion collisions (HIC) is 
strongly coupled [T] has provided motivation for under- 
standing the dynamics of strongly coupled non-Abelian 
plasmas through the gauge/string duality [2] (see [3] for 
a review of applications to the QGP). The simplest ex- 
ample of the duality is the equivalence between four- 
dimensional TV = 4 SU{N C ) super Yang-Mills (SYM) the- 
ory and IIB string theory on AdS$ x S 5 . Here we extend 
this example to the case in which the SYM plasma is spa- 
tially anisotropic. For accessibility by a broad audience, 
details will appear elsewhere [4]. Previous holographic 
studies of anisotropic plasmas include [5J [5] . One impor- 
tant difference with the gravity solution of [5] is that the 
latter possesses a naked singularity, whereas our solution 
is completely regular. 

Part of our motivation comes from the fact that the 
QGP created in HIC is anisotropic. An intrinsically 
anisotropic hydrodynamic description has been proposed 
to describe the early stage after the collision [7j. After 
that stage each little cube of QGP is isotropic in its own 
rest frame, but even in this phase certain observables may 
be sensitive to the physics in several adjacent cubes. 

Weakly coupled plasmas, both Abelian and non- 
Abelian, are known to suffer from instabilities in the pre- 
sence of anisotropies [5]. It is therefore interesting to 
understand whether this also happens in strongly cou- 
pled anisotropic plasmas. Our gravity solution exhibits 
instabilities reminiscent of weak-coupling instabilities. 

At a more theoretical level, motivation is provided by 
a connection with the fluid/gravity correspondence [5] 
and the blackfold approach to black hole dynamics [TUj . 
both of which assert that the effective theory describing 
the long- wavelength dynamics of a black hole horizon is 
a hydrodynamic theory. Inclusion of conserved p-form 
charges on the gravity side leads to anisotropic hydrody- 
namics pTj . 



The IIB supergravity solution that we will present is 
a finite-temperature generalization of that of |12j and: 
(i) it is static and anisotropic; (ii) it possesses a horizon 
and it is regular on and outside the horizon; (Hi) it obeys 
AdS§ x S 5 asymptotic boundary conditions. Staticity is 
required for simplicity, since e.g. we would like to study 
the thermodynamics of the system. The presence of a 
horizon is dual to the existence of a finite-temperature 
plasma in the gauge theory. Regularity guarantees that 
calculations are unambiguous and well defined. The 
boundary conditions ensure that holography is on its 
firmest footing and that the solution is solidly embed- 
ded in string theory. 

As in 12J, we deform the SYM theory by a ^-parameter 
that depends linearly on one of the three spatial coordi- 
nates, 9 — 2im- D7 z, where constant with dimen- 
sions of energy. In other words, we add to the SYM 
action a term 8S oc J 0(z) Tr F A F. The system we are 
describing is therefore a static plasma in thermal equi- 
librium in the presence of an anisotropic external source. 
Yet, translation invariance is preserved, since integration 
by parts yields SS oc -n D7 f dz A Tr (A A F + |A 3 ) . 

The dual gravity description is as follows. Since the 8- 
paramctcr is dual to the IIB axion \ |23j , we expect that 
in the gravity solution this will be of the form \ = az. 
It turns out [3] that a = An D7 /4irN c , where A = g^r M N c 
is the 't Hooft coupling. Since the axion is magnetically 
sourced by D7-branes, the solution can be interpreted 
in terms of a number of D7-branes wrapped on the S , 
extending along the xy-directions and distributed along 
the z-direction with density n D7 = dN or /dz [HIH]. For 
this reason we will refer to a and/or n D7 as a 'charge 
density'. Thus in the gravity description it is clear that 
isotropy is broken by the presence of anisotropic extended 
objects. Since their full backreaction is incorporated, the 
D7-branes are completely 'dissolved' in the geometry, just 
like the N c D3-branes that give rise to AdS$ x S 5 . Unlike 
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the case of flavour D7-branes [13], the D7-bran.es consid- 
ered here do not extend in the radial direction. Conse- 
quently, they do not reach the AdS boundary and they 
do not add new degrees of freedom to the SYM theory. 

As in |12j . the solution can be viewed as a renormali- 
zation group (RG) flow between an AdS geometry in the 
ultraviolet and a Lifshitz-like geometry in the infrared. 
At T = the Lifshitz metric (in string frame) possesses 
a naked curvature singularity [3], but this is hidden be- 
hind the horizon at T > 0. 

2. Solution. The ten-dimensional solution is a direct 
product, one of whose factors is an S 5 of constant radius 
L in the Einstein frame. Therefore it can be viewed as 
a solution of five-dimensional supergravity with cosmo- 
logical constant A = — 6/L 2 . Since only the metric g, 
the axion x, and the dilaton cf> are excited, it suffices to 
consider the axion-dilaton-gravity action 



(1) 

where we have set L = 1 and S GH is the Gibbons-Hawking 
boundary term. The Einstein-frame metric is 



,1s 2 = e -^- ( -TBdt 2 + dx 2 + dy 2 + Hdz 2 + ^ 



,(2) 



and x = az - Isotropy in the xy-directions is clearly re- 
spected, but not in the z-direction unless % = 1. The 
axion induces the anisotropy. J 7 is a 'blackening factor' 
that vanishes at the horizon, u — u n . The boundary is 
at u = 0. The dilaton only depends on the radial co- 
ordinate u, as do J 7 , B, and "H, which are completely 
determined in terms of <f>. This in turn obeys a third- 
order ordinary differential equation which we solved nu- 
merically [4]. The temperature is determined from the 
requirement that the Euclidean continuation of ^ be 
regular, and the entropy density from the area of the 
horizon. These quantities are well defined since the so- 
lution is static, i.e. the dual plasma is in thermal equi- 
librium (see Section 5). Fig. [ljleft) shows the entropy 
density as a function of a/T, normalized by the isotropic 
value s°(T) = tt 2 N 2 T 3 /2 pH]. This provides us with the 
following check. We see from the log-log plot that for 
small a/T the points lie on the horizontal axis, while for 
large a/T they lie along a line with slope 1/3. Thus at 
T » a we recover the isotropic result, whereas at T <C a 
we recover the Lifschitz scaling s oc a 

i/3 T 8/3 found 

m 

[12) . This interpolating behaviour is expected from the 
interpretation of the solution as an RG flow. 

3. Holographic stress tensor. The energy density 
and the pressures can be obtained from the holographic 
stress tensor, whose calculation requires the addition of 
counterterms to ([I]) . These can be obtained from P~5] [T5] 



and (in Euclidean signature) take the form 

S ^ = ^J d ' x Vl (z-\e 24, d lX d l x^-\ogv J d A x^yA, 

(3) 

where v is the Fefferman-Graham (FG) coordinate, 7 
is the induced metric on a v — vq surface, and the 
limit Vq — > is understood. A(^fij, <f>, \) is the conformal 
anomaly, which when evaluated on our solution takes the 
value A(^ij,(j),x) = N 2 a 4 /48ir 2 . 

From the results of [H] the stress tensor is seen to be 
diagonal, (Ty) = diag(_E, P±,P±,Pu), and to obey 



d i {T ij ) = 0, (T*) = A, 



(4) 



thus confirming that translation invariance is preserved. 
P± (P\\) are the pressures in the x, y (z) directions. As 
a consequence of the anomaly the transformation of the 
stress tensor under a rescaling of a, T contains an inho- 
mogeneous piece [HE], i-e. 



(Tij (ka, kT)) = k 4 (Tij (a, T)) + fc 4 log k A h 
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(5) 



where hij = diag (1,-1,-1,3). In turn, this means that 
the stress tensor must take the form 



(Tijia, T)> = a 4 ty {a/T) + log (a/ M ) A h 



(0) 



where the arbitrary reference scale fi is a remnant of 
the renormalization process, much like the subtraction 
point in Quantum Chromodynamics (QCD). Different 
choices of [i are simply different choices of renormali- 
zation scheme. We emphasize that the presence of this 
scale implies that the physics depends on the two dimen- 
sionless ratios T / \i and a/ \±, not just on T/a. Represen- 
tative plots of the energy and the pressures are shown in 
Fig. [ij center) . 

4. Thermodynamics. As usual, the free energy 
F(a,T) — E — Ts = —P± is obtained from the on- 
shell Euclidean action and satisfies (8F/dT) a = — s 
[H [TT]. Unlike the entropy density, which is scheme- 
independent, the energy density and the pressures are 
scheme-dependent (i.e. depend on (/,), but the thermody- 
namic relations among them are scheme-independent [3] • 
We recall that the necessary and sufficient conditions for 
local thermodynamic stability are 

c a = T(dS/dT) a >0, F"= (d 2 F/da 2 ) T >0. (7) 

5. Phase diagram. Approximate analytic solutions 
can be found in the limits T 3> a, fi and T <C a, /i, and 
these suffice to draw the qualitative phase diagram shown 
in Fig. fright), which we have also verified numerically 
[1]. F"(a,T) is negative in Zone I and positive in Zones 
II and III. P\\(a,T) - P°(T) is negative in Zones I and 
II and positive in Zone III, with P°(T) = ir 2 N 2 T 4 /8 the 
isotropic pressure. Note that each of the three zones in- 
cludes points with T — as well as points with arbitrarily 
large a and T. 
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FIG. 1: (Left) Entropy density as a function of a/T. (Center) Energy and pressures as functions of T/a for fixed a ~ 2.86 and 
\ogfj, = 1/2. (Right) Qualitative phase diagram. 



It follows that the homogeneous phase with uniform 
D7-brane density is in stable thermal equilibrium in Zone 
III. In particular, as can be seen from the positive slope 
of the continuous red curve in Fig. [TJ center), the specific 
heat is c a > everywhere. Also, the pressures and the 
energy are all monotonically increasing functions of T at 
fixed a, so the speed of sound in all directions is real and 
positive. There are no thermal instabilities anywhere in 
the phase diagram. 

In contrast, the homogeneous phase is in unstable ther- 
mal equilibrium against infinitesimal charge fluctuations 
in Zone I, where the second condition in ^ is violated. In 
Zone II the system is in metastable thermal equilibrium, 
since it is unstable only against finite charge fluctuations: 
the pressure in the z-direction is smaller than the pres- 
sure of the isotropic phase, and thus bubbles of isotropic 
phase can form and grow, forcing a compression of the 
charge in the z-direction [21] . In other words, in Zones 
I and II a carefully prepared homogeneous system with 
initial (a, T) will fall apart into a mixed phase consisting 
of high-density anisotropic 'droplets' or 'filaments' sur- 
rounded by isotropic regions [25] , The local charge den- 
sity a' > a will be the same in each of the droplets, and 
the pressure will exactly equal that of the isotropic phase 
at the same final temperature, P\\(a! ,T') = P°(T'). The 
pair (a',T') therefore lies on the continuous blue curve 
of the phase diagram. 

In Landau's theory of phase transitions the homoge- 
neous phase in Zones I, II and III would be described by 
a saddle point of the free energy with at least one unsta- 
ble direction, by a metastable local minimum, and by a 
stable global minimum, respectively. 
6. Discussion. Our system is in anisotropic thermal 
equilibrium. This is not surprising, since in the gauge 
theory isotropy is broken explicitly by an anisotropic 
external source. The string description makes it clear 
that the resulting system can be thought of as a fluid 
with a conserved, isotropy-beaking, two-brane charge 
(see e.g. pi]). 

It is remarkable that our solution is completely regular 
despite the fact that it incorporates the full backreaction 



of the D7-branes, whose number scales as n D7 ~ N c /X. 
Relatedly, we note that the parameter controlling their 
backreaction, Xn D7 /N cl is coupling-enhanced as in the 
case [18] of flavour D7-branes. 

The physics in Zones I and II shares some similarities 
with that of QCD at low T and finite baryon density 
[T§] . In that case the pressure of a chirally broken homo- 
geneous phase with density lower than a critical density 
no is negative (except in a tiny region of very small densi- 
ties). This indicates an instability towards the formation 
of 'droplets' of higher density no in which P = and 
chiral symmetry is restored, surrounded by empty space 
with n = and P = 0. In our case, the role of the chirally 
restored phase is played by the anisotropic phase, the 
analogue of n is a (see Fig. fright)), and the 'droplets' 
correspond to the regions of non-zero D7-brane density. 
These similarities suggest that the transition from the 
mixed phase to the homogeneous phase may occur via a 
percolation mechanism, as in some QCD models [20] of 
chiral symmetry restoration [26] . 

The instabilities we have uncovered are reminiscent of 
instabilities of anisotropic weakly coupled plasmas [5]. 
Somewhat pictorially, the main similarity is the tendency 
to 'filamentation', which in weakly coupled plasmas can 
be understood (very roughly) as the tendency of simi- 
larly oriented currents to cluster together. We emphasize 
though that there are obvious differences. In a weakly 
coupled plasma the anisotropy is 'dynamical' since it 
arises from the momentum distribution of the particles 
that compose the plasma. In contrast, in our case the 
plasma is static and intrinsically anisotropic because of 
the presence of dissolved extended objects. In any case, 
we stress that whether a real connection exists between 
the instabilities studied here and those of weakly coupled 
plasmas is a question beyond the scope of this letter. 

The instabilities of our solution are also reminiscent 
of the phenomenon of cavitation, i.e. the formation of 
bubbles of vapour in regions of a flowing liquid in which 
the pressure of the liquid drops below its vapour pres- 
sure. Cavitation has been proposed [H] (see also [2"2] ) 
as a mechanism that would lead to fragmentation of the 
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QGP into droplets that would subsequently evaporate, 
thus providing a new scenario for how hadronization is 
achieved. In that context the analogue of vapour pressure 
is the pressure of the vacuum, P = 0, whereas in ours it 
is the pressure of the isotropic phase. As above, how- 
ever, we emphasize that in the case of [H] the pressure 
drop is due to a dynamical effect, namely to the viscosity 
corrections that result from the expansion of the plasma. 
In contrast, in our case this is a static effect presumably 
resulting from the interaction of the extended objects in 
the plasma. 

Note that an instability discovered in [T^] does not 
directly apply here, since Lifshitz (instead of AdS) boun- 
dary conditions were assumed in |12j . 
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